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We consider a network of topological defects which can partly decay into neutrinos, photons, 
flj ■ baryons, or Cold Dark Matter. We find that the degree-scale amplitude of the cosmic microwave 

[j, ' background (CMB) anisotropies as well as the shape of the matter power spectrum can be con- 

siderably modified when such a decay is taken into account. We conclude that present predictions 
concerning structure formation by defects might be unreliable. 
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I. INTRODUCTION 



Two challenging paradigms to explain structure formation in the Universe are currently developped, namely cos- 

mological inflation ||l| and topological defects [g|. On the one hand, inflation is a simple theory, based on the linear 

rS^ evolution of acausal and coherent initial perturbations produced during an accelerated phase of expansion of the early 

M-i Universe. Several public codes are available, notably CMBFAST |g|, to compute the CMB anisotropies of a given 

Q model in few minutes, and these predictions are robust. On the other hand, topological defects, which are supposed 

5-H to have formed after a cosmological phase transition in the early Universe, are much more difficult to handle, because 

c/2 , of the highly non- linear structure of their dynamics. After some pioneering work by Bennett, Bouchet, Stebbins [H, 

^ ' Shellard, Allen Q, Perivolaropoulos |^, Caldwell et al. 0|, the first detailed predictions of some defect models have 

^ [ been published only recently. In particular, Turok, then Pen, Seljak and Turok [p| have considered global defects, 

JJ"! ■ as well as Durrer and collaborators M who studied more carefully the large TV limit. Following Vincent, Hindmarsh 

rN \ and Sakellariadou ||I3|, Battye, Albrecht and Robinson 11^, have studied a network of line- like segments with given 

S ■ correlation length and velocity distribution, as well as Pogosian and Vachaspati ||l^ who also considered the wiggly 

" " ' structure of the strings. Allen, Shellard and collaborators |^^, and Contaldi, Hindmarsh and Magueijo |lj] have 

considered local cosmic strings. 

Several ground experiments (Saskatoon, Python V, TOCO, see [|l^; see also ||l^ for an up-to-date compilation of 
the current results) have by now probed the degree-scale anisotropy of the microwave sky. They seem to indicate the 
presence of a high peak in the spectrum. Now, in contradiction with these observations, topological defect models 
do not produce generically much more power on the degree scale than on the largest scales observed by the COBE 
satellite [Q. The reason is that, although the defects evolve according to causal processes, thus producing power on 
small angular scales on the last scattering surface, one must also take into account their gravitational interaction with 
the photons after last scattering. This so-called Integrated Sachs- Wolfe effect (see e.g. |1^), known to be negligible in 
most inflationary scenarios, greatly increases in the case of defects the power on large angular scales, thus contradicting 
the present observations (see however |1^, where the wiggly structure of strings is properly taken into account, yielding 
more power on the degree scale). Moreover, in most defect theories, the situation is also worsened by the fact that 
scalar, vector and tensor modes all contribute significantly to the overall CMB anisotropies, whereas only the scalar 
contribution is expected to produce more power on small angular scales than on larger ones. 

It must however be stressed that, in all the numerical defect models already mentionned, an important physical 
effect, to wit their decay into gravitational radiation and/or elementary particles, has been considered in most cases 
in a phenomenological way only (but see however [|l9|). Turok et al. and Durrer et al. have imposed and checked 



that their global defect stress-energy tensor is conserved. In cosmic string numerical simulations one has to deal with 
the problem of loop production and evolution, which has been treated in various ways. Turok et al. have introduced 
an extra radiation fluid into which their cosmic string loops are supposed to decay, Battye et al. have introduced 
an extra fluid with given constraints on its stress-energy tensor, and Magueijo et al., whose stress-energy tensor is 
not conserved due to the fact that the smallest loops are extracted out of the simulation, "dump" the string energy 
losses into cither extra-fluids with no anisotropic stress and various equations of state, or in one or the other of the 
existing background fluids, that is the baryons, photons, neutrinos or CDM. Finally, Shellard et al. treat the loops as 
relativistic point masses, and Pogosian et al. treat them as small segments. 

The aim of this paper is first to show that the microphysics of defects imposes that energy be released in the cosmic 
fluids and that the predictions concerning the CMB anisotropics change drastically when even a small fraction of the 
network energy is released directly in the form of photons or baryons, and, to a lesser extent, neutrinos. As for the 
matter power spectrum, we will see that it is greatly affected if the network energy is released into photons, neutrinos, 
or baryons. 

Such a lack of generic predictive power for seeds-based models has already been pointed out in another con- 
text by Durrer and Sakellariadou |20|| , who have used various ansatze for the stress-energy tensor components of 
scaling coherent seeds and have found a wide variety of results for the positions and heights of the induced CMB 
anisotropics Doppler peaks. One could argue that these results were obtained by using very ad-hoc hypothesis for 
the seeds correlators, but at least three other works have reached similar conclusions by using more realistic models. 
Perivolaropoulos |g] has pointed out that a correct modelling of the gravitational interaction between cosmic strings 
wiggles and electrons could induce more small scale anisotropics than previously expected, and this analysis was con- 
firmed by the more precise numerical works of Pogosian and Vachaspati [|2| . Finally, the inclusion of decay products 
in the specific model of Contaldi et al. |Q proved to have an important effect on the CMB anisotropics spectrum. 
Interestingly enough, these two models which add some small-scale physics give a better agreement with observations 
than those which do not include any small-scale physics. We show here that such effects are not restricted to any 
specific model but hold for a very large class of defects and decay processes. 

We shall first discuss the microphysics behind such effects, then show how one can generically take them into 
account, and finally present some results. 

II. THE DECAY PRODUCTS OF REALISTIC STRING MODELS 

Topological defects fall in various classes |g]. 

Uncharged Goto-Nambu strings form loops which are usually supposed to decay dominantly into gravitational 
radiation |2|,Q , thus guaranteeing the scaling behaviour of the network (that is the fact that the ratio of the energy 
density contained in the network with the total energy density is constant in time) . This extra gravitational radiation 
added to the background is observationally constrained by e.g. the millisecond pulsar timing measurements p3| . 

Uncharged global strings, when seen from a distance, appear very much like local ones, but with an energy per 
unit length renormalized to include long-range interaction effects pj]. The overaU global string network is therefore 
expected to behave much like a Goto-Nambu network, except for the small loops |g5[ . The main difference between a 
local and a global string loop is that while the former is supposed to radiate mostly in gravitational waves, the latter 
looses energy mainly through massless Goldstone boson radiation, a process known to be far more efficient ||25| . This 
new extra component added to the background is only constrained by nucleosynthesis pffl. 

The last potentially interesting class of strings in the context considered here is that of Grand-Unified (GUT) 
superconducting cosmic strings [ p7| (only GUT strings can be relevant in large scale structure formation and CMB 
fluctuations). Again, their network evolution is almost the same as that of Goto-Nambu strings because most of the 
evolution takes place when the currents flowing along the strings are negligible pq |. This is however true only for 
infinite strings or for those having very large radius of curvature. 

When a superconducting loop decays, its energy can be very efficiently released in the background directly in the 
form of electromagnetic radiation ^^ . The electromagnetic radiation cannot propagate because of the surrounding 
plasma, the influence of which cannot be neglected. The waves form shells where energy is concentrated. These shells 
are the basis for the explosive Ostriker, Thompson and Witten large scale structure formation model pfl] and yield 
a distorsion in the microwave background by modifying its spectrum, implying non zero values for the parameters /i 
(chemical potential distortion, see e.g. pll) and y (characterizing the Sunyaev-Zel'dovitch effect, see e.g. 132]). This 
very stringent constraint, together with e.g. nucleosynthesis constraints [ p3| almost rules the scenario out, leaving as 
the only possibility that the radiation can only be emitted at much higher frequencies (to allow propagation in the 
plasma). This requires that the strings can be considered current-free until they have shrunk to a sufficiently small 
size, the wavelength of the radiation being proportional to the emitting loop radius. As a result, one can consider the 



network evolution as, again, that of a Goto-Nambu network with the difference that part of the energy contained in 
the loops can now be transferred directly into the photon fluid. 

When such a loop shrinks however, the integrated current being conserved, the energy it contains per unit length 
increases and its effects become more and more important on the string dynamics. The resulting loop distribution can 
accumulate to stationary states known as vortons [p3 . This overproduction of vortons breaks the scaling behaviour 
of the network, thus leading to such a cosmological catastrophe that one must assume that the vortons are sufficiently 
unstable to decay in less than one Hubble time. 

The new phenomenon one therefore needs to consider is the fate of the small charged loops. As discussed above, 
they can decay into gravitational radiation, or directly into photons. The third, largely overlooked possibility when it 
comes to compute CMB anisotropics or matter power spectrum, is that they might decay into their constituents [g5| , 
namely massive Higgs and gauge bosons and the particles they couple to and that make the current, collectively 
referred to as "X" -particles. As these particles have masses comparable to the grand unification scale, their decay 
products can only be estimated by the relevant QCD extrapolations at high energy |3^. The standard scenario is 
that the primary X-particle will decay into a lepton (usually an electron) and a quark which subsequently initiate 
a hadronic shower. Once the shower has evolved, one ends up with roughly 3% nucleons [p6| , and the rest in tt's, 
which, because of their decays or interactions with other background fluids, turn into neutrinos, photons and electrons. 
These interactions have even been used to try and explain the ultra high energy cosmic ray ||37| enigma by means 
of topological defects |38| . Note that this picture does not take into account the possibility that part of the decay 
product be a stable particle, i.e. a constituent of the dark matter. The Lightest Super Particle is a possible candidate, 
for instance if SUSY is demanded, as should be the case for any GUT model. 

This analysis of small superconducting loops has led us to conclude that they mostly decay (hence ensuring that 
their network scales) directly into the constituents of the universe rather than into extra-fluids such as gravitational 
radiation. 

In fact one can also argue that non superconductiong strings, because they intercommute, can also partly decay 
into photons, baryons, neutrinos or dark matter, and not only into gravitational radiation or Goldstone bosons for 
the global ones. Even global defects can also decay |g] : when the gradients are strong enough high energy particles 
and gravitational radiation are produced. One might expect however that the energy losses in the form of background 
fluids be far more efficient for superconducting, rather than uncharged, defects. 

III. A TWO-PARAMETERS DECAY MODEL 

When one solves the classical (Goto-Nambu) equations of motion for the string network in a Friedmann-Robertson- 
Walker (FRW) background, completed by a set of rules which fixes the intercommutation of long strings into loops, in 
the stiff approximation fS^] framework (in which the defect network energy is a first order perturbation to the FRW 
background) , and if no extra physics is added to the problem, the corresponding stress-energy tensor of the strings 
must be covariantly conserved. If it is not, that only means that the numerical integration is not precise enough. The 
same holds true for global defects. 

Now, as we have seen in § |ll| , extra physics must be added. Indeed intercommutation leads to the formation of more 
and more small loops per Hubble volume which would prevent scaling to take place and soon dominate the evolution 
of the universe, were they not eliminated by turning into gravitational radiation and/or various elementary particles. 
The string stress-energy tensor cannot then be conserved. This is particularly clear in numerical simulations when 
loops smaller than a given size are phenomenologically extracted by hand from the network (as in e.g. [|l4| ). 

In the semi-analytic approach we use here (see below) , we model active seeds by a stress-energy tensor O^^, which, 
from the start, embodies the scaling properties of the network : first, deep in the radiation-dominated or matter- 
dominated eras, its equal-time correlators (0^j,(x, ?7)0^^(x', 77)) depend on (x — x')/r] only, second they vanish at 
small scales (77 is conformal time and a;* are comoving coordinates). Therefore this effective stress-energy tensor 
describes in fact, not only the defects themselves but also their decay products, as long as these decay products are 
themselves "seeds", that is an active perturbation added to the background fiuids. Now among the possible decay 
products described in §||, only the gravitational and Goldstone radiations enter in that category. Hence, if the defects 
were scaling only via the production of this "extra" radiation, the effective stress-energy tensor we use, being the 
sum of the stress-energy tensor of the defects and that of the extra radiation they produce, would be conserved. In 
contrast, in numerical simulations where the defect stress-energy tensor is not conserved, such decay products must 
be treated as an extra fiuid with given equation of state. 

However, as we have seen, in the case of superconducting strings at least, loops decay mainly through the emission 
of high energy particles which soon turn into (mostly) neutrinos and photons, which must be added to those already 
existing in the background. Therefore we shall describe the defects by an effective stress energy-tensor which will not 



be conserved in order to take into account the decay of the loops into high energy particles, and we shall dump the 
string energy losses not in an extra component but into the background fluids. 

In the semi-analytic approach initiated by Durrer and collaborators Q, the defect network and the extra "seeds" 
they decay into are modelled by a stress-energy tensor O^^ which acts as a source for the linearized Einstein equations 
and induces inhomogcncitics in the cosmic fluids. The ten components of this stress-energy tensor are then drastically 
constrained by a number of physical requirements |40|-|42|] , which take into account that the seeds : 

1. are statistically homogeneous and isotropic, 

2. are created at a phase transition in an up to then perfectly homogeneous and isotropic universe (hence obeys 
specific causality pffl and matching mU conditions), 

3. evolve deep in the radiation era and deep in the matter era in a way which is statistically independent of time 
(scaling requirement p2|). 

One can now suppose that the defect network is made of long strings of comoving curvature radius Rl > 1/Ti. 
{TL is the comoving Hubble parameter), which eventually interconnect so as to form loops. When they become too 
small these loops decay at a certain rate (depending on the details of the decay processes). The loss of energy of the 
stress-energy tensor can therefore, in a rough approximation, be modelled as 

i?,,er = F\ (1) 

F"" = -fiQ"/ with fi = TY{k - kc) (2) 

where D^ is a covariant derivative (we assume fiat spatial section), a hat denotes a Fourier transform, Y is the 
Heavyside function, k oc 1/ Rl is the comoving wavenumber. F is a free function which determines the decay rate of 
the loops into the decay products, and kc is the scale above which the loops begin to decay. This F'^ is then injected 
into the other fiuid perturbation equations with branching ratios Xf , so as to ensure that the total stress-energy tensor 
is conserved (see appendix H for the actual equations we solve) . The branching ratios Xf are in principle calculable 
when the decay microphysics is explicited. 

IV. RESULTS 

A. CMB anisotropies 

We have computed numerically the CMB anisotropies and the matter power spectrum such decaying seeds produce. 
We first summarize our findings before going into discussing them. Dumping energy into CDM has negligible effect 
on the CMB anisotropies spectrum. Dumping energy into neutrinos has a relatively small influence on the spectrum, 
which is modifled by a few tens of percent. This can either increase or decrease the amplitude of the spectrum [we hav e 
found that this depends on the parameters of the decay term (g), as well as on the ansatze for the correlators (Aq,A7)]. 



On the contrary, the spectrum is tremendously affected by injecting energy into either photons or baryons, although 
the effect is much stronger when one injects energy directly into photons, which can boost the degree-scale amplitude 
of the spectrum by a factor as large as 50. The precise amplitude of the boost depends of course on the underlying 
defect model, but all those we have checked present the same qualitative features. 

A simple interpretation of these results is that energy, when injected into either neutrinos or CDM, influences the 
photons perturbations only through gravity, so that the influence of this energy injection on the CMB anisotropies is 
small. On the opposite, by dumping energy into baryons or photons, one directly affects the evolution of the perturbed 
photons density and/or velocity since photons and baryons are strongly coupled through Thomson diffusion before 
recombination. The angular scale at which this effect is strongest is given by the angular size of the "decay scale" 
(i.e. the Hubble radius when kc = Ti) at the last scattering surface, that is ~ 1 deg. 

For example in Fig. |l|, we consider the scalar contribution to the CMB anisotropies of a specific cohere nt model 



where we have chosen the equal time correlators for the pressure and the anisotropic stress [see eqns. (|A6|,A7)] in such 
a way that they respect the standard causality and scaling requirements (see appendix for more details). We have also 
chosen the relative amplitudes of the correlators so that, when the stress-energy tensor is conserved, we "mimic" the 
situation of more realistic incoherent numerical defect models (i.e. the CMB anisotropies spectrum lacks power on the 
degree-scale with respect to observations; the coherent model presented here exhibits acoustic oscillations which are 
expected to be smoothed out by decoherence, see ^). The two parameters of the decay model, kc and F [see eq. (||)] 
have both been chosen equal to Ti. We have also studied many other models such as the "pressure model" Q] [where 



ones sets the anisotropic stress 11'* to 0, see eqns. (A6 A7)], or the "anisotropic stress model" (where the pressure f"* 
is set to 0, see e.g. [Q), which aU yield similar qualitative results. In all what follows, the CMB data points are taken 
from |1^ and the matter power spectrum data points are those of the APM catalog |Q and of Dekel et al. |4^ . 

In Figs, y and ^ we study the influence of the parameters kc and F for the model of Fig. |l| where energy is equally 
released into photons or neutrinos. By increasing kc (Fig. 0), which gives essentially the wavenumber scale at which 
energy is injected, we get that, as expected, the spectrum is enhanced on smaller angular scales, i.e. the bump is 
shifted to the right. On the contrary, by increasing F (Fig. H), namely the defects decay rate, energy is transferred 
more rapidly into the background fluids, and therefore at larger wavelengths. The spectrum is thus enhanced on 
larger angular scales, so that the bump is shifted to the left. 

B. Matter power spectrum 

Dunrping energy into CDM has also negligible effect on the (baryonic) matter power spectrum. On the contrary, the 
spectrum is strongly affected by injecting energy into either photons or neutrinos. It has the consequence of reducing 
the excess of energy on small scales, because of the free streaming of these relativistic particles (see e.g. Q). Finally, 
injecting energy into baryons gives an intermediate result. Fig. ^ summarizes these results by showing the matter 
power spectrum corresponding to the model of Fig. ^. 

C. Observational constraints 

In all these numerical analysis, we have assumed that the radiation emitted by the high energy decay product is 
immediately thermalized as soon as it is produced. However this is only an approximation. Three unwanted effects 
can be caused by this radiation. First, the radiation emitted before the last scattering surface may not have had 
the time to thermalize, thus leading to a distortion of the CMB black-body spectrum; this would be the case if too 
much energy is injected between z ~ 10^ and z ~ 10'^. However we inject a small amount of energy compared to that 
already contained in the CMB (roughly 10~^ per Hubble time since one is in the radiation-dominated era), so that 
(considering the current observational bounds on the fj, and y parameters) one can assume it had time to thermalize, 
see pM. Second, the radiation emitted after the last scattering surface will not be thermalized, and therefore could 



produce an important 7-ray background |33|. Third, the high-energy particles could photo-dissociate the He nuclei. 



thus producing lighter nuclei such as He and D, which in turn can produce too much Li |48| . 

These three effects are already constrained by the observations of the CMB spectrum, the diffuse 7-ray background 
and the light elements abundances, but our scenario still happens to be tenable, as shown on Fig. ^ (taken from |p3|). 
In the opposite case, this would disproof most scenario of structure formation seeded by topological defects which 
dominantly decay into photons. 

V. CONCLUSION 

We have included some microphysics, up-to-now largely overlooked, in the description of topological defects. This 
microphysics deals with the decay of defects, and notably superconducting cosmic strings, into background fluids, 
rather than gravitational radiation as usually assumed. This decay has important observational consequences, which 
may (depending, of course, on the exact interaction one considers) put the defect models in a better position when 
confronted to the current observational data. 

One could of course argue that the simple model we have considered here (coherent seeds, crude interaction term) 
is too naive, however our purpose was mainly to illustrate the consequences of this idea rather than to study a more 
specific realistic model, which we keep for later work. Since the main consequences of dumping energy into the 
background fluids do not seem to depend too much on the details of the equal time correlators we have considered, 
but rather on the details of the interactions, we strongly advocate for a more careful study of the branching ratios. 

Let us recall that even in inflationary models the inclusion of microphysics and interaction between fluids is ab- 
solutely crucial in order to make accurate predictions. For example, if one neglects these by not solving the exact 
Boltzmann equation for the photons and/or the neutrinos, or by not solving the accurate kinetic recombination 
equations, one finds an large excess of power at small angular scales (see e.g. Fig. 4 of Q). 

It is therefore clear that predictions of topological defect models concerning structure formation not only require 
today's state-of-the-art heavy detailed numerical simulations, but also a rigorous description of their non gravitational 
interactions before reliable conclusions can be drawn. 
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APPENDIX A: CMB ANISOTROPIES CALCULATIONS 

We assume standard cosmological parameters : flat Universe without cosmological constant, baryon density and 
Hubble parameter such that fib = 0.05, h — 0.5, three massless neutrinos and standard recombination. 

The stress-energy tensor of the seeds can be decomposed as a sum of scalar, vector and tensor |B3| random fields. 
The scalar part, the only one we shall consider here, can be written in terms of four random fields p*, u", P^ and tt^, 
as : 

KO^oo ^ /5^ (Al) 

KQ'o, = -ihv\ (A2) 

K% = P'6,j + Qfc2% - k,k^ ff . (A3) 

{K is Einstein's constant.) To compute the CMB anisotropics, we need the two-point correlators of the stress-energy 
tensor. The homogeneity of the distribution imposes : 

e^,(k,r?)e^,(k',77')) =5(k-k')C',,.p.(k,r;,V), (A4) 

where the k-dependence of C is fixed by the requirement that the distribution is isotropic. The correlators can be 
decomposed as sums of "coherent eigenmodes" |8| : 



C',..p.(k,ry,,7') = J] A«c«(k,,7)c«(k,,/), 



(A5) 



(») 



where the CpllCk, rj) are the correlators of coherent sources. The four random fields that describe coherent sources are 
all proportional to the same normalised random variable e(k) and hence are described by four random fields. Two 
of these functions are constrained by eqns. (p]-0), and the two other must be either imposed by hand or determined 
using some more detailed modelling. We take simple ansatze for the pressure P'* and the scalar anisotropic stress fl'* , 
in practice (up to a normalization constant) : 



U'' = -4:r]-^k'^r]^) exp{-k^r]^)e{k). 



(A6) 

(A7) 



This (rather arbitrary) choice of correlators satisfies the requirements imposed by causality, scaling, etc 142], and is 
tailored in order to obtain a CMB anisotropics spectrum which lacks power on the degree scale. Of course, the precise 
shape of these functions can in principle be obtained by numerical simulations, and the CMB anisotropics are known 
to depend very much on the shape of these free functions |2^, but, as already stressed, we have carefully checked that 
the effect we are interested in does not qualitatively depend on them. 

We solve (in the flat-slicing gauge) the well-known (see e.g. M,B3,b2l) linearised Einstein equations which couple 
the seeds network to the cosmic fluid inhomogeneities, corrected by the decay term (yj^). They read : 






7 



5c — k^Vc + Xcfl 



Si, = k'^vij + Xb/i 



37^2^c 
P' 



«7 = -4^7 



3 



/if 



— k fl-y — 4* — $ — K {vj — -Of,) -I- Xjfi 



AH^n^ 



k = -\s. + lk^ti.-^-^ + ^.P^^, 

Vc ~ —HVc — $ -f Xcfl 



iW^Vlc ' 



~ 4 p v^ 

Vb = -Uvb - $ - -— /i {vb - v^) + ^fcA o-i-/2o ' 

6 pb SH'^ilb 



(A8) 
(A9) 

(AlO) 
(All) 
(A12) 



where fi/, 5/, Vf, 11/ are respectively the density parameter, the density contrast, the velocity and the anisotropic 
stress pertubations for the fluid /, the subscripts 6, c, v, 7 mean respectively the baryonic. Cold Dark Matter, neutrino 
and photon fluid, ^ and $ are the two Bardeen gravitational potentials, k is the photons differential opacity [pTl , and 
a dot denotes a derivation with respect to the conformal time. The terms in /i account for the seeds decay, and the 
branching ratios Xf are constants such that : 



/ 



1. 



In addition, the two equations for the Bardeen potentials are : 



m - Ti")] * = 37^2 ^ ^f [6f - 3(1 + iOf)vf 
f 

!> = * - 37i;2 ^ fifiVftif ~ ti". 



p' -3Hv', 



Finally, the photons and neutrinos fluids obey a Boltzmann equation which reads (for i > 2, see e.g. 152] or 
more details) : 



^ 2^+1 



£Ai-' -{£ + l)A^+i ) -k(A'- SiM , 



^ 



21 
1 



- (Vl - 4/i?2^^-i -{1+ 1)^1-4/(^+1)2^^+1) - K (P^ + ^P-^h,^ 



^^=iol^ 



^ Vei;^ 



A^ = 



2i + 



-(M:-i-(£+l)Af/i), 



(A13) 

(A14) 

(A15) 

g for 

(A16) 
(A17) 
(A18) 
(A19) 



with Ai and Ei being respectively the £-th moment of the temperature and electric-type polarisation distribution 
functions of the species /, and where Pf is the coupling between temperature and polarization. The lowest multipoles 
of the distribution function are related to the perturbed stress-energy tensor components by : 

(A20) 
(A21) 

(A22) 



Sf-- 


.4A0, 


Vf-- 


-I^J 


n/ = 


^S-? 



Eqns. (Al- A22 ) are solved numerically using a Boltzmann code developped by one of us (A.R.) using standard initial 
conditions p:l|. The CMB anisotropics at present conformal time tjq are then calculated by the linc-of-sight integration 
method : 



Al(fc,7?o)= / ke 



Vo 



-d^ + ^ + <P]]i{k{T]a-i^))d7] 



ke "^ {-kvb) j'^{k{T]Q - ri))dr] 



Vo 



* + $ j^(fc(77o-?7))dr7 



1)0 5 

ke-'^-P^ (3j;'(fc(77o - v))+Je{k{vo - v))) dry. 



(A23) 



In the linc-of-sight integration formula ( A23| ), we have omitted the terms arising from the energy injection due to 
the seeds decay into photons. These terms are proportional to e^^/i and therefore are generated between the last 
scattering surface and today. The photons emitted because of this decay have therefore very high energy and do not 
have time to thermalize. Hence, they do not participate to the microwave background, but rather to the diffuse 7-ray 



background, which might put some interesting constraints on these models as discussed in §IVC 



Finally, the scalar part of the 2-point correlator of the CMB anisotropics is decomposed as 

^^(2^+l)CfP,(cos0), 



ST 6T 

^(ni)^(n2) 



ni .n2— cos(0) 



(A24) 



where P£ is the Legendre polynomial of order £, and the coefficients C/ are deduced from the photon distribution 
multipoles by : 



C 



|Ai(fc,r,o)|^fc^dfc 



(A25) 
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FIG. 1. CMB anisotropies in a model where the defect stress-energy tensor is not conserved, with kc = T = Ti.. The solid 
line represents the case with conserved stress-energy tensor, as well as the case where one dumps energy into CDM, which 
are almost identical. The short-dashed and long-dashed lines represent the cases where one dumps energy into neutrinos and 
baryons respectively. The highest (dotted) line represents the case where energy is injected into photons. The dot-dashed line 
(which fits the data points best) shows that by tuning by hand the branching ratios (here, x^ = x-, = 0.5 and Xc ~ xt — 0) it 
is possible to be in much better agreement with the data points. Note that the precise values of the best-fit branching ratios 
depend on the model one considers. 
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FIG. 2. Influence of the decay scale kc on the CMB anisotropies. We use the "best fit" model of Fig. hi, with 50% of the 
energy relased into photons and another 50% into neutrinos, and we vary the parameter kc from Ti. to 57i. The solid line 
represents as in Fig. hi the case where the seed stress-energy tensor is conserved, and the case kc = T — Ti corresponds to the 
dot-dashed curve, also plotted in Fig. hi 
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FIG. 3. Influence of the decay rate F on the CMB anisotropics. We use the "best-fit" model of Fig. hi, with 50% of the energy 
relased into photons and another 50% into neutrinos, and we vary the parameter F from Ti. to 5Ti.. The solid line represents 
as in Fig. hi the case where the seed stress-energy tensor is conserved, and the case kc = T = H corresponds to the dot-dashed 
curve, also plotted in Fig. hi. 
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FIG. 4. Matter power spectrum in the same model as in Fig. y, with kc — T = H. The solid line represents the case 
where the seed stress-energy tensor is conserved. Dumping energy into neutrinos (short-dashed line) or photons (dotted line) 
has almost the same effect, whereas dumping energy into CDM (solid line) has again negligible effect. Dumping energy into 
baryons (long-dashed line) boosts the spectrum on larger scales, whereas it damps it on smaller scales. Note that the sharp 
bump around k — 10""^ hMpc~^ shows the scale at which one injects energy (here, the Hubble radius since kc = Ti). 
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FIG. 5. Maximal energy release in units of the CMB energy density allowed by the constraints from the observed diffuse 
7— ray background at 5 GeV (dotted curve), CMB distortions (dashed curve), and ''He photo-disintegration (solid curve) as a 
function of redshift z (taken from Sigl et at). These bounds apply for instantaneous energy release at the specified redshift 
epoch. The lowest (dot-dashed) curve shows the energy release in our model assuming that ~10% of the energy is released 
into photons (energy releases into the other fluid do not add to these constraints). Higher energy release into photons might 
violate the 7-ray background measurements. Note that the strongest constraint arise today from the diffuse 7-ray background, 
but better measurements of the CMB spectrum could give an even stronger constraint. 
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